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Abstract 

We compute the crossover exponent <j> describing the crossover from the 
random-exchange to the random-field critical behavior in Ising systems. For 
this purpose, we consider the field-theoretical approach based on the replica 
method, and perform a six-loop calculation in the framework of a fixed- 
dimension expansion. The crossover from random-exchange to random-field 
critical behavior has been observed in dilute anisotropic antiferromagnets, 
such as Fe a .Zni_a;F2 and Mn r Zni_. r F2, when applying an external magnetic 
field. Our result <j) = 1.42(2) for the crossover exponent is in good agreement 
with the available experimental estimates. 

PACS Numbers: 75.10.Nr, 75.10.Hk, 64.60.Ak 



1 



Typeset using REVTeX 



The most studied experimental realizations of random-field Ising systems are dilute 
anisotropic antiferromagnets in a uniform magnetic field applied along the spin ordering 
axis [1-3]. A simple lattice model is provided by the Hamiltonian 

H = J ^2 PiPj SiSj - Pi s h (1) 

<ij> i 

where J > 0, the first sum extends over all nearest-neighbor sites, Sj = ±1 are the spin 
variables, and pi are uncorrelated quenched random variables, which are equal to one with 
probability p (the spin concentration) and zero with probability l—p (the impurity concen- 
tration). In the absence of an external field, i.e. H = 0, the critical behavior of dilute Ising 
systems (above the percolation point of the spins) belongs to the universality class of the 
random-exchange Ising model (REIM), which differs from the standard Ising one. See, e.g., 
Refs. [4,5] for recent reviews. The applied uniform field H gives rise to a different critical 
behavior corresponding to the universality class of the random-field Ising model (RFIM). 
See e.g. Refs. [6-8] for recent reviews on RFIM. 

For t = (T — T N )/T N — > 0, T N being the Neel temperature in zero field, and H — > the 
singular part of the free energy can be written as [1] 

F=\u t \ 2 -<*f(H 2 \u t \-*), (2) 

where u t ~ t + a\H 2 + a^t 2 is the scaling field associated with the temperature, a is the 
specific-heat exponent of the REIM, f(x) is a scaling function, and is the corresponding 
crossover exponent. As a consequence of the crossover scaling (2), the critical temperature 
in the presence of the external field H is given by 

T C {H) ^T N + cH 2 ^ - aH 2 , (3) 

for sufficiently small H. Experimental measurements on dilute antiferromagnets yielded 
rather precise estimates: = 1.42(3) obtained for Fe :c Zni_ :c F2 [9], = 1.43(3) for 
Mn^Zni^Fs [10], and = 1.41(5) for Fe^Mgi^Cls [11]. 

On the theoretical side, a computation of was presented in Ref. [12] using the field- 
theoretical e expansion. We recall that in the REIM the e expansion is actually an expansion 
in powers of \fe. The exponent was computed to 0(e) [12], obtaining 

0/ 7 = 1 + 0.16823e 1/2 - 0.22666e + 0(e 3/2 ), (4) 

where 7 is the REIM susceptibility exponent. Apart from showing that 0^7, where 
7 1.34 (see, e.g., Refs. [13-15]), this series does not yield a reliable estimate of for three- 
dimensional systems. This is generically true for the y/e expansion of any critical exponent 
of the REIM [4,5]. 

In this paper we determine the crossover exponent using an alternative field-theoretical 
method based on a fixed-dimension expansion in powers of appropriate zero-momentum 
quartic couplings, which we compute to six loops. As we shall see, the analysis of the series 
provides a quite precise estimate of 0, 

0=1.42(2), (5) 
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in good agreement with experiments. This result was already anticipated in Ref. [16]. 

The field-theoretical approach is based on an effective Landau-Ginzburg- Wilson (/? 4 
Hamiltonian that can be obtained by using the replica method [17], i.e. 



r l N l N 

/ d M o E [i^if + rtf] + T ]T (u + voSij) tftf}, (6) 
J z i=i ^- ?j=i 



where </?j is an A^-component field. The critical behavior of the REIM is expected to be 
described by the stable fixed point of the Hamiltonian in the limit iV — > for u < [18]. 
The most precise field-theoretical results for the critical exponents have been obtained by 
analyzing the fixed-dimension expansion in powers of the zero-momentum quartic couplings 
u, v related to uq,vq (i.e. u = Z u uo/m and v = Z v vo/m where Z UjV = 1 + 0(u,v)), which 
have been computed to six loops [15,19]. We refer the reader to Ref. [15] for notations and 
definitions. 

In the presence of a spatially uncorrelated random field h(x) with zero average and 
variance h 2 , one may still apply the replica method. Averaging the replicated random field 
term 

J d 3 xh(x)J2<Pi(x) (7) 

i 

over the Gaussian distribution of the field, one obtains a new term proportional to 

h 2 I ^xY,Vi{x)Vj{x), (8) 

that must be added to the Hamiltonian 7i v 4 [12]. This term causes the crossover from the 
REIM to the RFIM. Thus, the REIM-to-RFIM crossover exponent is determined by the 
— > limit of the renormalization-group (RG) dimension y T of the quadratic operator [20] 

Ti3=tpMj, i^j. (9) 

In order to evaluate yx and the corresponding crossover exponent <f) = yr^, we define a 
related RG function from the one-particle irreducible two-point function with an 
insertion of the operator Tjj, i.e. 

^r\^)ij,ki — Zj, 1 Aij k i, (10) 
where A^ k \ = 5ikSji + Su5jk so that Z T (0) = 1. Then, we compute the RG function 



, , d\nZ T 

Vt{u,v) = — 

a mm 



n d\nZ T d\nZ T 

Uo,VQ 



where f3 u and /3 V are the /5-functions. The exponent i]t is given by the value of t)t(u,v) at 
u = u*, v = v*, where (u*,v*) is the REIM stable fixed point. Finally, the REIM-to-RFIM 
crossover exponent is obtained by using the RG scaling relation 

= (2 +f] T - 7]) v = 7 + T] T U, (12) 
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where 7, v, and r] are REIM critical exponents. 

We computed the function r^ T,2 ^(0) to six loops. The calculation is rather cumbersome, 
since it requires the evaluation of 563 Feynman diagrams. We handled it with a symbolic 
manipulation program, which generates the diagrams and computes the symmetry and group 
factors of each of them. We used the numerical results compiled in Ref. [21] for the integrals 
associated with each diagram. The resulting six-loop series of t)t(u,v) is 

rj T {u,v) = -\u + ±u 2 + ±uv - 0.0357673 m 3 - 0.0483240 u 2 v - 0.00421548 uv 2 (13) 
+0.0343748 u 4 + 0.0762616 u 3 v + 0.0416943 u 2 v 2 + 0.00896116 uv 3 
-0.0408958m 5 - 0. 121377 u A v - 0.104778m 3 v 2 - 0.0373712 u 2 v 3 - 0.00527015 uv A 
+0.0597048 u 6 + 0.227662 u 5 v + 0.287108 u A v 2 + 0.17231 u 3 v 3 
+0.0552124 u 2 v 4 + 0.00759654 uv 5 + 

where u and v are the rescaled couplings u = u/(6ir) and v = 3v/ (167r), and the dots indicate 
higher-order terms. In our analysis we also considered the series <p(u, v) corresponding to 
0, which can be obtained by using Eq. (12) and the series of the RG functions ^(u, v) and 
v{u,v) corresponding to the critical exponents 7 and z/, cf. Ref. [15]. 

In order to obtain an estimate of 7]t and 0, the corresponding six-loop series must be 
resummed and then evaluated at the fixed-point values u*, v*. Although the perturbative 
expansion is not Borel summable [22] , various resummation schemes have been proposed and 
employed, obtaining rather reliable results for the critical exponents, see, e.g., Refs. [4,5]. 
We analyzed t]t(u,v) and (j)(u,v) by using the resummation methods outlined in Ref. [15] 
and the estimates u* = —18.6(3) and v* = 43.3(2) obtained by Monte Carlo simulations 
in Ref. [13], which turn out to be more precise than those obtained from the zeroes of the 
/3-functions computed to six loops. 

We obtained t]t = 0.095(30) from the analysis of i]t(u,v), = 1.43(3) from <p(u,v), and 
4> = 1.42(2) from l/(f)(u,v). The errors are related to the spread of the results obtained 
by using different resummation methods and different resummation parameters within each 
method, see Ref. [15] for details. Using the available estimates of the critical exponents 
(e.g., 7 = 1.342(6), v = 0.683(3) and r] = 0.035(2) from Monte Carlo simulations [13,14], 
7 = 1.330(17), v = 0.678(10) and r\ = 0.030(3) from the analysis of the six-loop field- 
theoretical expansions [15]), the above- reported result for r/ T and the RG scaling relation 
(12) give = 1.41(2). From these results we arrive at the final estimate reported in Eq. (5), 
which is in good agreement with the available experimental results obtained for various 
uniaxial anisotropic antiferromagnets. 
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